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Abstract 

The correlation function of two pseudoscalar charmed quark currents 
with a positive hadronic spectral density is employed to obtain upper 
bounds on the decay constants of D and D a mesons. Including all known 
terms of the operator-product-expansion of this correlation function in 
QCD and taking into account the estimated uncertainties, we obtain jo < 
230 MeV and fo, < 270 MeV. Comparison with the decay constants 
determined from D — > Ivi and D a — > Ivi measurements, reveals a tension 
between the bound and current experimental value of /d s . 



1 Introduction 



The decay constants of charmed D + and D s mesons, defined via the hadronic 
matrix elements: 

(0|d7 M 7 5 c|L» + (p)) = ifopn , (0|s7 Al 7 5 c|L» s (p)} = if Ds P» , (1) 

are extracted from the branching fractions of purely leptonic decays D + — > l + i>i 
and D s — > l + vi (I — \i, r), respectively. Recent CLEO measurements of these 
decays yield [TJ [2] : 

/z, = 205.8 ± 8.5 ± 2.5 MeV, (2) 
f Da = 267.9 ±8.2 ±3.9 MeV, (3) 

assuming \V c d\ — \V US \ and |V CS | = \V u d\- Other experimental results for fjj and 
fn s are consistent with the above intervals (see [3] for a review). Note that the 
observed SU (3) /(-violation in D s and D decay constants turns out to be larger 
than in the light pseudoscalar mesons where fx = 155.5 MeV and /„• = 130.4 
MeV @]. 

The decay constants of heavy mesons are accessible in lattice QCD. The 
recent results: f D = 207 ± 4 MeV and f Da = 241 ± 3 MeV [5], obtained with 
the number of sea-quark flavours Nf = 3 have quite small errors. The lattice 
value of fr> is in a good agreement with $2$, whereas fn s is smaller than ([3]). 
This puzzling situation caused discussions of a possible non-standard physics 
(see e.g., [SJ[7]). The previous N f = 3 lattice QCD result 0: f D = 201 ± 3 ± 17 
MeV, fn s — 249 ± 3 ± 16 MeV, as well as the two recent Nf — 2 calculations: 
f D = 205 ± 7 ± 7 MeV, f Ds = 248 ± 3 ± 8 MeV [9] and f Ds = 257 ± 3 ± 3 ± 5 
MeV (preliminary) [10] quote larger errors, still their central values for fo s are 
smaller than in ([3]). Recent results for /d w in quenched lattice QCD can be 
found in [TT1 [TJ] and a discussion of the accuracy of lattice determinations in 

An alternative way to calculate /d (s) is provided by QCD sum rules |T3 . The 
method is based on the operator-product expansion (OPE) of the two-point cor- 
relation function evaluated in deep spacelike region, in terms of perturbatively 
calculable Wilson coefficients and QCD vacuum condensates. The result of 
OPE is related with the hadronic matrix element (JT|) via dispersion relation and 
quark-hadron duality. The QCD sum rule calculation of the heavy-meson (-D( s ) 
and -B(s)) decay constants has a long history, starting from [Ml [15j [16] . More 
recent calculations using QCD sum rules are in [TTl [18l [19] and finite-energy sum 
rules in [50] ■ A review of earlier results can be found in [3T]. Not going into 
further details, let us only mention that QCD sum rules yield fp a > fp, and 
the predicted intervals for both decay constants are in the ballpark of the lattice 
QCD results, albeit with larger uncertainties. In particular, a "systematic" un- 
certainty of QCD sum rules which is difficult to quantify in a model-independent 
way, is caused by the quark-hadron duality approximation. Having in mind a 
possible confrontation with experiment, an update of the QCD sum rule pre- 
dictions for /d (s) is a timely task which is however beyond our scope here. 

The aim of this paper is to remind that the same correlation function which 
is used to obtain the QCD sum rule, provides upper bounds on fo and fn s ■ 
These bounds simply follow from the positivity of the hadronic spectral density 
of the correlation function and are independent of the quark-hadron duality 
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approximation. In what follows, we will calculate the upper bounds for D and 
D a meson decay constants and find that the bound on fjj a is on the verge of 
disagreement with the experimental value. 



2 Derivation of the bounds 

We start from the correlation function of two charmed pseudoscalar quark cur- 
rents j5 = (m c + mdjdi^c, the divergence of the axial-vector current in (fT]): 



u( q 2 ) = i { dW<o | r{i,(«)jt(o)} | o) = £ (o\m!MM . (4) 

h=D,D*Tr,. h y 



For definiteness, the channel of the charged D meson is considered, the cor- 
responding correlation function Tl s (q 2 ) for D s -channel is obtained by a simple 
replacement of the light-quark flavour d — » s. On the r.h.s. of the corre- 
lation function, via unitarity condition, is written as a sum over all hadronic 
states with D quantum numbers, a very schematic representation of the disper- 
sion integral over hadronic spectral density. The ground _D-state contribution 
to ([3]) contains the square of the decay constant: (0|j5 = font^. 

At large virtualities, q 2 <C m% the correlation function ^ is dominated 
by short distances. In that region we approximate n(g 2 ) by OPE, including 
the contributions of the perturbative quark-loop diagrams (with gluon radiative 
corrections) and the terms with vacuum condensates. The latter are suppressed 
by inverse powers of m 2 c — q 2 and it is sufficient to include condensates up to 
dimension d < 6. For virtual c and light quarks propagating in the correlation 
function, it is convenient to adopt the MS scheme. The most complete expres- 
sions for OPE in this scheme are presented in [T7] where they were used for 
the QCD sum rule calculation of /b and fs s ■ The corresponding expressions 
for -meson case are simply obtained by replacing the heavy quark mass 
m& — > m c , hence, there is no need to represent here the explicit formulae. We 
only write down a schematic form of the OPE result after Borel transformation: 



n(M 2 )= I rfs (v)" P (n) ^ s/M2 

"=0,1,2 , / 

+ E (v)" n Si (M2)+ £ Hd(M2) - (5) 

n=0,l d=4,5,6 

The Borel parameter M 2 replaces q 2 and represents the virtuality scale at which 
the correlation function is calculated. In the above, the perturbative contribu- 
tions of quark-loop diagrams are represented in a convenient form of dispersion 
integrals over spectral densities, p(s) = ilmll(s). The a s -expansion of p(s) in- 
cludes the two-loop [H] and three-loop [H] gluon radiative corrections. For the 
latter the program 'Rvs.m', kindly provided by the authors of [52] was used. The 
second line in (0 includes the contributions of the quark condensate with the 
0(a s ) correction calculated in [17], as well as the standard gluon-, quark-gluon- 
and 4-quark-condensate terms of OPE denoted by the respective condensate 
dimension d = 4, 5, 6. 
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Borcl-transforming ([¥]), we use the positivity of the hadronic spectral den- 
sity. Hence, H( S \(M ) provides an upper limit for the ground-state .D( s )-meson 
contribution to the hadronic sum in (U . The resulting upper bound for the _D/ S ) 
decay constant is: 

m% /M 2 

n (s) (M2) e J' . (6) 

where H^(M 2 ) is calculated from The bound is valid at any M 2 , at which 
one can trust the OPE ([5|) . Naturally, one has to find the most restrictive value 
of the r.h.s. in ©. 

3 Numerical analysis 

Numerical results for the bounds are obtained adopting the c-quark mass in- 
terval fh c (fh c ) — 1.29 ± 0.03 GeV which covers the recent determinations [23] 
from charmonium sum rules with 0(a^) accuracy, and, conservatively, we dou- 
ble the uncertainty. Note a good agreement of this determination with the 
recent lattice QCD result [21] for fh c . For the strange quark mass we adopt 
m s (2GeV) = 98 ± 16 MeV, the interval of QCD sum rule determinations with 
0{a A s ) accuracy [25J. Fixing m s and using the ChPT relations [2BJ we obtain 
m„(2GeV) = 2.8 ± 0.6 MeV and m d (2GeV) = 5.2 ± 0.9 MeV. The quark-gluon 
coupling and quark masses are taken with 4-loop running, employing the pro- 
gram provided in [27], with a s (mz) = 0.1176 ± 0.002 [J. We assume equal 
renormalization scales for the quark masses and a s . The quark condensate den- 
sity is obtained using Gell-Mann - Oakes - Renner relation: (0|gg|0)(l GeV) = 
-mlfl/[2{m u +m d ){\ GeV)] = -(250±^ MeV) 3 , (with U = (130.4±0.04±0.2) 
MeV and = 139.57 MeV g]). 

The ratio of the strange and nonstrange quark condensate densities (0|ss|0) = 
(0.8 ± 0.3)(0|gg|0), as well as the intervals for the quark-gluon, gluon and four- 
quark condensate densities not quoted here for brevity, are taken as in |29j . 
where a different correlation function using these universal parameters was cal- 
culated (see also [30] for a review of vacuum condensates). The suppression 
of d = 4, 5, 6 terms in OPE makes uncertainties of the respective condensate 
densities inessential for our numerical results. 

The upper bounds on fjj and fn s calculated from © are plotted in Fig.l 
as a function of M 2 at central values of all input parameters. The default 
renormalization scale is \i = 1.5 GeV, close to /i ~ \J m 2 D — m 2 used in sum rule 
calculations. In MS scheme the perturbative expansion works reasonably well: 
0{a s ) and 0(a 2 ) terms in ([5]) are, respectively, < 30% and < 10% of the total 
perturbative contribution. 

As expected, the bound for Jd s is larger than the one for fr>. Both bounds 
grow and become less restrictive at larger M 2 , implicitly due to increase of the 
relative weight of the continuum and excited states in the Borel-transformed 
hadronic sum. The most restrictive upper bounds 

f D < 220 MeV, f Ds < 250 MeV , (7) 

1 Positivity bounds of various form for this correlation function can be found in the literature 
starting from |28| . 
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Figure 1: The upper bounds for frj (left) and for fu^ (right) as functions of the 
Borel parameter squared at \i = l.5GeV (solid) and [i — 3.0 GeV (dashed). 



are reached around M 2 = 1.2 GeV 2 . At these values of the Borel scale, the 
sum of power suppressed d — 4, 5, 6 condensate contributions is less than 4% 
of the total H(M 2 ), hence the condensate expansion can be trusted. The scale 
dependence turns out to be mild, e.g., increasing the scale from fi = 1.5 GeV 
to fi = 3 GeV, one obtains the most restrictive bounds at M 2 ~ 1.4 GeV 2 , and 
they are only slightly (by about +5 MeV) shifted (see Fig.l). Decreasing the 
scale up to fi = 1.0 GeV, produces a more pronounced shift, by about -15 MeV, 
so that the bounds are even more restrictive. However, the NNLO correction 
at this scale reaches ~ 20% of the total perturbative contribution, signalling 
that the perturbative expansion is less convergent numerically. Hence, to be on 
a conservative side, we use the bounds obtained at the default scale, including 
their uncertainty (at fixed M 2 ), caused by the variation 1 GeV < /i < 3 GeV, 
in the error budget discussed below. 

The bounds have uncertainties caused by the limited accuracy of QCD pa- 
rameters (quark masses, a s and condensate densities) and by the scale-dependence 
The individual uncertainties are estimated by varying all inputs one-by-one 
within the adopted intervals. The results are collected in the table: 



Variation of input 


/u-bound uncertainty 


fo 3 -bound uncertainty 


m c 


± 2.0% 


±2.0% 


m s 




± 1.4% 


Us 


±0.7% 


± 0.7% 


<o|??|o> 


-2.5 /o 


-1.9 /o 


(0|ss|0}/<0|g<?|0) 




±5% 


d=4,5,6 condensates 


± 1.0% 


± 1.0% 


scale \i 


-2.1 /o 


-3.0 /o 


total in quadr. 


±4.8% 


±6.9% 



Summing up separate uncertainties in quadratures, we obtain ~ ±10 MeV 
(±20 MeV) for the f D (f Ds ) bound in J7J). The uncertainty of the /^-bound 
is naturally larger, due to the variation of m s and the spread in the ratio of 
strange and nonstrange condensates. Conservatively, we shift the bounds (J7|) 
up by their respective total uncertainties, yielding our final estimate: 



f D < 230 MeV, f Ds < 270 MeV . (8) 
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4 Discussion 



The lattice results [5J [HI GO EH EQ HI] for /d and /d s are consistent with the 
upper bounds ©■ The decay constant of D-meson ^ extracted from experi- 
ment also obeys its bound. Formally, the experimental value Q of fr> s does not 
contradict ([S]) as well, especially if one takes into account the experimental er- 
rors. However, it seems quite unnatural for this bound to be almost completely 
saturated by the D s contribution. In the correlation function Tl s (q 2 ) at timelike 
q 2 there are intermediate hadronic states located above D s , starting from D*K 
{D*tt is forbidden by isospin symmetry), and including the radial excitations 
of D s and continuum hadronic states with D s quantum numbers. It is difficult 
to evaluate in a model-independent way the size of individual contributions of 
these states to the hadronic spectral density. Their Borel-transformed sum, a 
positive quantity, can be estimated, subtracting the £> s -contribution from the 
correlation function H S (M 2 ) calculated from OPE: 

£ (° I h\h)(h\j$ I 0)e"^/ M2 = H s (Af 2 ) - flm^e-^J^ , (9) 

h=D*K,... 

and substituting the experimentally determined fn 3 from ([3]). For normalization 
we divide both parts of this equation by the D s contribution and calculate the 
ratio: 

_E (0\jB\h)(h\jt\0)e-</ M2 
R s (M 2 )= h = D * K >- (10) 

fl) s m D s e ° s 

from ([9|), as well as the analogous ratio R(M 2 ) for the D-meson channel where, 
correspondingly, the experimental result for fn is used. Note that in the SU(3) fi 
symmetry limit the correlation functions for D s and D channels and their 
hadronic components are equal: H S (M 2 ) = H(M 2 ), fr> s = /d, mD = ™_d,, 
R S (M 2 ) = R(M 2 ), etc. Our calculation of IT S (M 2 ) and II(A// 2 ) explicitly takes 
the 5 , C/(3)/;-violation effects into account, via differences of s- and u,d- quark 
masses and qq (q — u, d) and ss condensates. Since here we also use quite 
different values of fo^ and fo, the results for the hadronic quantities R S (M 2 ) 
and R(M 2 ), are drastically different, as seen from Fig. 2. At 1 < M 2 < 2 
GeV 2 the hadronic sum for D s is even negative, violating unitarity. At larger 
A/ 2 , the share of higher states in the Z? s -channel correlation function is strongly 
suppressed with respect to the same characteristics for Z?-channel. Maximizing 
the uncertainties in OPE (as we did it in obtaining the conservative bounds ([8])), 
one can shift both curves upwards, making R s marginally consistent with posi- 
tivity, however the large SU (3) ^-difference remains. A natural question is then: 
why is SU (3) ^-symmetry so strongly violated in this correlation function? 

In general, the reliability of the bounds ^ depends on the convergence of 
OPE. Including the 0(a 2 ) perturbative loops and condensate contributions up 
to d — 6 we obtain a reasonable numerical convergence. One, however, cannot 
completely exclude, that some nonperturbative, e.g., instanton-like effects con- 
tribute to II(Af 2 ) beyond OPE, and modify the bounds. However, even if such 
effects influence the convergence of OPE, it is hardly possible that they simulta- 
neously produce a drastic SU(3) /^-violation. Moreover, a similar OPE in terms 
of perturbative loops and power-suppressed local condensates for the light-quark 
(pseudo)scalar currents was successfully used to obtain the lower bounds for the 
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Figure 2: The lower (upper) solid curve is the sum of higher-state hadronic 
contributions to the Borel-trans formed correlation function in D s (D ) channel, 
normalized by the D s (D ) contribution, as a function of the Borel parameter 
squared. The OPE for the correlation function II S (M 2 ) (H(M 2 )) is calculated 
at central values of input parameters and the experimental value for /d s (fr>) 
is used. Dashed lines indicate the experimental errors in the decay constant 
determination. 

light-quark masses [3T], up to 0(af ) accuracy [3^. If there were noticeable ef- 
fects beyond OPE in the correlation function of charmed pseudoscalar currents, 
their role would have become more pronounced for the light-quark mass bounds, 
obtained with a light quark replacing the virtual c quark. 

5 Conclusion 

Concluding, we find that the current experimental result for fo 3 only marginally 
obeys the upper bound for this decay constant obtained from the OPE for the 
two-point correlation function. This is in contrast to the D-meson case where 
the analogous bound is consistent with fo inferred from experiment. A more 
precise knowledge of c- and s-quark masses and of the ratio of strange and non- 
strange quark condensates will make the bounds more accurate. Furthermore, 
the contributions of higher states to the correlation function in D s channel esti- 
mated using the experimental value of fu s , are strongly suppressed as compared 
with the corresponding contributions for the D meson channel. 

If the noticed difference between strange and nonstrange charmed meson 
channels remains in future, from QCD point of view, that could indicate a 
presence of unaccounted nonperturbative effects in the pseudoscalar correla- 
tion function of charmed quarks which violate not only OPE, but also SU(3)fi 
symmetry. However such effects seem not to manifest themselves in the lattice 
calculations. 
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